We report on high resolution spectra of the bright QSO 
HS 0747+4259 (Zem= 1-90, V = 15.8) observed to search 
for intermediate redshift O VI absorption systems. The spec- 
tra were obtained by means of the Space Telescope Imaging 
Spectrograph (STIS) at the Hubble Space Telescope (HST) and 
the High Resolution Echelle Spectrometer (HIRES) at the W. 
M. Keck telescope. We identify 16 O VI systems in the range 

1.07 < z < 1.87. Among them, six systems with Zabs = 1.46 - 

1 .8 exhibit a sufficient number of lines of different ionic transi- 
tions to estimate the shape of the ionizing radiation field in the 
range 1 Ryd < i? < 10 Ryd. All recovered UV ionizing spec- 
tra are characterized by the enhanced intensity at i? > 3 Ryd 
compared to the model spectrum of Haardt & Madau (1996). 

Jjf^ This is in line with the observational evidence of a deficiency of 
strong Lya absorbers with A^(H l) > 10^^ cm~^ at z < 2. The 
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UV background shows significant local variations: the spec- 



tral shape estimated at z — 1.59 differs from that obtained at 
Q z = 1.81 and 1.73. A possible cause of these variations is the 
presence of a QSO/AGN at z ~ 1.54 — 1.59 close to the Une of 
sight. No features favoring the input of stellar radiation to the 
ionizing background are detected, limiting the escape fraction 
of the galactic UV photons to /esc < 0.05. 
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Abstract. 



1. Introduction 

The ionization state of the intergalactic medium (IGM) is main- 
tained by the metagalactic ionizing radiation field. The current 
paradigm considers the metagalactic UV field to be produced 
by the integrated radiation of QSOs reprocessed by the clumpy 
IGM (Haardt & Madau 1996, hereafter HM96; Fardal et al. 
1998). Model calculations of the spectral energy distribution 
(SED) of the UV background depend crucially on the distribu- 
tion of the Lya forest lines - the number of absorption lines per 
unit redshift and per unit H I column density. This distribution 
is stiU poorly known. For instance, there is observational evi- 
dence suggesting a deficiency of strong absorbers with A^(H i) 
> 10^'' cm^^ at z < 2 compared to higher redshifts (e.g. 
Kim et al. 2001). These absorbers are the main source of the 
He II continuum opacity and their deficiency may result in a 
ionizing UV background much harder at energies E > 4 Ryd 
than model predictions based on biased absorber statistics (see 
HM96, Sect.5.14). Another fact is that at 2: < 2 we observe sig- 
nificant spatial variations in the absorption line number density 
(Kim et al. 2002). 

Closely related to the SED of the background ionizing radi- 
ation is the problem of a possible galactic (stellar) contribution. 
It is supposed that the emissivity of QSOs decreases dramati- 
cally at z < 2, whereas the input of galactic radiation to the 
ionizing background becomes more significant (e.g., Bianchi et 
al. 2001; Bolton et al. 2005). This conclusion is made from es- 
timates of the mean intensity of the ionizing background at 912 
A. However, a putative galactic contribution would affect not 
only the intensity at the hydrogen ionization threshold but also 
the shape of the ionizing spectrum making it softer at A < 304 

Send ojfprint requests to: S. A. Levshakov, 
lev @ astro.ioffe.rssi.ru 

* Based on observations obtained at the W. M. Keck Observatory 
and the HST 



A (Leitherer et al. 1999; Giroux & ShuU 1997). Thus, recover- 
ing the SED in the range E > A Ryd will enable us to estimate 
the relative contributions of QSOs and galaxies to the ultravio- 
let background. 

To probe the shape of the ionizing radiation, optically thin 
absorption systems (10^^ cm"^ < A''(Hl) < 10^"^ cm~^) 
with metal lines can be used (Chaffee et al. 1986; Bergeron & 
Stasinska 1986). The ionization thresholds of ions usually ob- 
served in the intervening absorbers (Si Il-Si IV, C II-C IV, N III, 
N V, O VI) lie between 1 and 10 Ryd and, thus, these bound- 
aries limit the accessible energy range. Usually the shape of the 
ionizing continuum is estimated by a trial-and-error method. 
Recently we developed a special procedure based on tech- 
niques from the theory of experimental design that allows us 
to restore the SED directly from the analysis of absorbers with 
many metal lines (Reimers et al. 2005; Agafonova et al. 2005). 

In this paper we report on the SEDs obtained from the anal- 
ysis of metal absorbers with z^hs = 1. 46- 1.81 identified in the 
spectrum of HS 0747H-4259. All absorbers studied reveal lines 
of O VI. This ion can be produced either coUisionally in high 
temperature gas (Tkin > 10^ K) or by photoionization due to 
the background ionizing radiation. Particular interest in coUi- 
sionally ionized O VI stems from the fact that O Vi can trace the 
so-called 'warm-hot intergalactic medium' (WHIM) predicted 
in cosmological simulations (Cen & Ostriker 1999; Dave et al. 
2001). Thus, the inferred shape of the ionizing continuum may 
help us to understand which mechanism is preferable for the 
interpretation of the observed O vi. 

We also estimate the frequency of O vi absorbers in units 
of absorption distance, ANqyi/ ^X. 

This paper is organized as follows. The observations are 
described in Sect. 2. Sect. 3 contains a short description of the 
computational methods used to analyze the absorption systems, 
which are studied in detail in Sect. 4. The results obtained are 
discussed and sunmiarized in Sect. 5. 
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Table 1. Metal absorption-line systems toward HS 0747+4259 
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Mg 112796,2803 ; Fe II26OO 
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Note: O vi AA1031, 1037 lines from systems with 2abs > 1.8617 fall in the gap between the HST/STIS and Keck/HIRES spectra 



2. Observations 

The UV spectrum of the quasar HS 0747+4259 (0em= 1-90, 
V = 15.8; Reimers et al. 1998) was obtained with HST/STIS 
in September 2001 with a total exposure time of 54 180 seconds 
over 19 orbits. The spectrum was taken with the medium reso- 
lution NUV echelle (E230M) and 0.2" x 0.2" aperture, provid- 
ing a resolution of A/AA 30 000 (FWHM 10 km s"! ). 
The spectrum covers the range from 2136 A to 2970 A. The 
data reduction was performed by the HST pipeline. After addi- 
tional inter-order background correction, the 12 individual ex- 
posures were co-added. The resulting signal-to-noise S/N is 
3-10 per pixel, lower than expected because the quasar was 
fainter by a factor of 2 in the UV compared to lUE observa- 
tions (Reimers et al. 1998). 

Another two integrations of HS 0747+4259 were obtained 
with HIRES on the Keck-I telescope; 3600 seconds on 2001 
February 28, and 5400 seconds on 2001 March 01. The CCD 
was the single Tektronics engineering grade with 2048 by 2048 
24 micron pixels that was fixed in HIRES from 1994 to mid- 
2004. Both integrations used a single setup covering 3060 A 
to approximately 4560 A with no gaps between the 39 spec- 
tral orders. We used the C5 dekker which contains a 1.148 
arcsec wide slit hat gives high throughput and approximately 
8 km s~^ FWHM resolution. Before each integration of the 
QSO, 300 second integrations of the flux standard star G191 
B2B were obtained. A Thorium-Argon lamp spectrum was 
used for the wavelength calibrations. The spectrum was ex- 
tracted and reduced using Tom Barlow's MAKEE package, 
following procedures described in Kirkman et al. (2003) and 





cm 

NUI 


i,i 

NV 


Si I II 


OIII OIV OVl 

J \ 






' \ 

; \ 


■ MF (AGN) 

- HM96, z=1.8 




\ \ ^ 


UV background SI 




\\ \ 


UV background S2 

, 1 




\ \ 

V ■ 

\ : 
\ ■ 



1 10 

E, Ryd 



Fig. 2. Schematic picture of the z = 2 metagalactic (smooth line) 
and AGN-type (long-dashed line) ionizing backgrounds from Haardt 
& Madau (1996) and Mathews & Ferland (1987), respectively. The 
spectra are normalized so that Juihf = 1 Ryd) = 1. The positions 
of ionization thresholds of different ions are indicated by tick marks. 
The results of the restoring procedure are shown by the dotted (SI) 
and short-dashed lines (S2). For details, see text. 

Suzuki et al. (2003). The obtained S/N decUnes systematically 
into the UV from ~ 30 at 4500 A to ~ 5 below 3200 A. 
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Fig. 1. Hydrogen and metal absorption lines associated with the Zabs = 1.8073 system toward HS 0747+4259. Continuum normaUzed intensities 
are shown by dots, along with 1 cr error bars. The zero radial velocity is fixed atz = 1.8073. The smooth curves are synthetic spectra calculated 
with the ionizing background SI (Fig. 2). The spectra are convolved with the corresponding spectrograph line-spread function. The physical 
parameters are listed in Table 2, Col.3, 5, and 6 (subsystems A, B, and C, respectively). The dashed lines show synthetic spectra computed with 
the HM96 ionizing spectrum (Fig. 2), the physical parameters are given in Cols.2 and 4, Table 2. Bold horizontal lines mark pixels included in 
the optimization procedure. The synthetic profiles of unmarked absorption lines were calculated in a second round using the velocity v{x) and 
gas density nn (x) distributions already obtained in the optimization procedure. 



3. Computational methods 

The absorption systems are analyzed by means of the 
Monte Carlo Inversion (MCI) procedure described in detail 
in Levshakov, Agafonova & Kegel (2000, hereafter LAK), 
and with modifications described in Levshakov et al. (2002, 
2003a,b). Here we briefly outUne the basics needed to under- 
stand the results presented below in Sect. 4. 



The MCI is based on the assumption that all lines observed 
in the absorption system are formed in a continuous medium 
where the gas density, nH^x), and velocity, v{x), fluctuate from 
point to point (here x is the space coordinate along the line of 
sight). 

Further assumptions are that within the absorber the metal 
abundances are constant, the gas is optically thin for the ion- 
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izing UV radiation, and the gas is in thermal and ionization 
equilibrium. This means that the fractional ionizations of all 
ions are determined exclusively by the gas density (or, equiva- 
lently, by the ionization parameter U oc l/nn) and vary from 
point to point along the sightUne. Since the ionization curves 
(i.e., the dependence of the ion fraction on U) are different for 
different ions, the ionic line profiles are different (non-similar) 
as well. Another important fact is that for a given point within 
the line profile the observed intensity results from a mixture of 
different ionization states due to irregular random shifts of the 
local absorption coefficient (see Fig. 1 in LAK). 

The fractional ionizations are determined by the SED of the 
background ionizing radiation which is treated as an external 
parameter. 

The radial velocity v{x) and gas density nYi{x) are con- 
sidered as two continuous random functions which are repre- 
sented by their sampled values at equally spaced intervals Aa;. 
The computational procedure uses adaptive simulated anneal- 
ing. The fractional ionizations of all elements included in the 
analysis are computed at every space coordinate x with the pho- 
toionization code CLOUDY (Feriand 1997). 

In the MCI procedure the following physical parameters 
are directly estimated: the mean ionization parameter Uq, the 
total hydrogen column density A''h, the line-of-sight velocity, 
cTv, and density, cry, dispersions of the bulk material [y = 
nn{x)/nQ], and the chemical abundances Z^, of all elements 
involved in the analysis. With these parameters we can further 
calculate the column densities N^, for different species, and the 
mean kinetic temperature Ttin- The mean gas number density 
is related to the parameters of the gas cloud as [see eq.(28) in 
LAK] 

no = ^(l + aJ). (1) 

Here riph is the number density of photons with energies above 
IRyd 

n^h. = —r J912 / — — , (2) 

where c, /i, Vc, and are the speed of light, the Planck con- 
stant, the frequency of the hydrogen Lyman edge, and the spe- 
cific intensity (in ergs cm~^ s~^ sr~^ Hz^^). If the mean gas 
number density no is known then the line-of-sight thickness L 
of the absorber can be evaluated as L = iVn/no. 

An important issue in the analysis of the absorption systems 
is the treatment of unidentified blends which can affect the line 
profiles from the Lya forest. Since the method supposes that 
all ions trace the same underlying gas density and velocity dis- 
tributions, it is possible to reconstruct both distributions using 
the unblended parts of available Unes. To clarify which parts 
may be blended several test runs with different arrangements 
of lines are carried out until a self-consistent fit for all lines 
observed in the system is found. 

As mentioned above, the spectral shape of the ionizing radi- 
ation is treated as an external parameter: some standard ioniz- 
ing spectrum is selected, corresponding ion fractions are cal- 
culated and then the MCI analysis is carried out with these 
fractions inserted. If the selected spectrum fails to reproduce 



the observed line intensities or some other physical inconsis- 
tencies arise (e.g. odd element abundance ratios) the search 
for a more appropriate spectrum can be performed. The cor- 
responding computational technique is described in Reimers 
et al. (2005) and Agafonova et al. (2005). It is based on the 
response surface methodology from the theory of experimen- 
tal design and includes (i) the parameterization of the spectral 
shape by means of a set of variables (called 'factors'), (m) the 
choice of a quantitative measure (called the 'response') to eval- 
uate the fitness of trial spectral shapes, and {in) the estimation 
of a direction (in the factor space) which leads to a spectrum 
with better fitness. The optimal (best fitness) spectral shape is 
that which allows us to reproduce the observed intensities of 
all Unes detected in the absorption system without any physical 
inconsistencies. 

How well the spectral shape can be recovered depends on 
the number of metal lines involved in the analysis: the more 
lines of different ionic transitions of different elements are de- 
tected in the absorption system, the more constrained is the al- 
lowable set of shapes. When only a few lines are detected, they 
can be used to estimate the shape in some restricted energy in- 
terval. For example, lines of C iii, C iv and O vi allow us in 
some cases to estimate the depth of the He 11 break {E > A 
Ryd). Any other a priori information concerning absorbers un- 
der study should be considered as well in order to distinguish 
between possible solutions. 

4. Absorption systems towards HS 0747+4259 

In the spectrum of HS 0747H-4259 we identified over 50 ab- 
sorption systems, among them 25 with metal lines (Table 1). 

From these systems, only a few turned out to be suitable for 
analysis by MCI. The selection criteria included the presence 
of both lines in the C IV and O VI doublets, and at least one ad- 
ditional carbon line (e.g. C III) in order to fix the ionization pa- 
rameter. As an initial guess for the UV background, the HM96 
spectra at each corresponding redshift was used. The mean in- 
tensities of the ionizing background at A = 912 A, J912, needed 
to calculate the number density of photons and the line-of-sight 
size of the absorbing gas, were also taken from HM96. 

Note that these J912 values only slightly differ (within 
20%) from those obtained from the proximity effect at 2: = 
1-2 (Scott et al. 2002) and from measurements of the mean 
radiation flux transmitted through the Lya forest (Tytler et al. 
2004; Jena et al. 2005). 

All computations below were performed with laboratory 
wavelengths and oscillator strengths taken from Morton (2003) 
for A > 912 A and from Verner et al. (1994) for A < 912 A. 
Solar abundances were taken from Asplund et al. (2004). 

The uncertainties of the physical parameter estimates {Uq, 
Nil, o-y, (Ty, and Z^) are about 20%. 

4.1. Absorption system atz = 1.8073 

This complex absorption system containing lines of different 
ionic transitions is shown in Fig. 1. Unfortunately, some of 
these lines are corrupted by bad background subtraction and/or 
high noise (Oiii 832, Oiv 788, probably also Cm 977). Our 
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Table 2. Physical parameters of the z^hs = 1.8073 metal absorber towards HS 0747+4259 derived by the MCI procedure with the Haardt & 
Madau and modified UV background spectra (marked, respectively, by HM96 and SI in Fig. 2). Column densities are given in cm~^ 
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"uncertainty caimot be estimated since line is corrupted by background subtraction; 

'^column densities for putative O VI absorption at —70 km s~^ < v < —20 km s"^ and at 

30 km s~^ < w < 50 km s~^ are not given because of probable blending with some Lya forest Unes 

(observed intensities of O VIA1032 and A1037 are inconsistent); ''can be blended 



first computational runs revealed that the system consists of 
three subsystems with different physical parameters: low ion- 
ization subsystems at —90 km s"^ < w < km s~^ (subsystem 
A) and km s~^ < t; < 80 km s~^ (subsystem B) seen in 
C II-C IV, Si Il-Si IV and O II-O IV and overlapping high ion- 
ization subsystem(s) seen in CiV and Ovi at —20 km s^^ < 
V <20 km s~^ (subsystem C) and perhaps in O iv, O vi at 30 
km s~^ < u < 50 km s~^ and in O vi at -70 km s~^ < i; < 
-20kms-^ 

In the following the subsystems A, B and C are treated sep- 
arately. Calculations of subsystems A and B were carried out 
with the ionizing spectrum of HM96 at 2; = 1.8 (solid line in 
Fig. 2). The obtained physical parameters are given in Table 2, 
C0I.2 and 4. Due to low S/N it was possible to reproduce al- 
most all line profiles to within the noise level, i.e. the spec- 
trum of HM96 may be considered as consistent with the ioniza- 
tion state of both subsystems. However, some overestimation of 
Si IV 1 393 in the subsystem A and underestimation of C 11 1334 
and Si 11 1260 in the subsystem B hint to probable inadequacy 
of the adopted spectral shape (Fig. 1, dashed lines). Low S/N 



and bad background subtraction do not allow us to restore the 
SED of ionizing radiation via the directed search procedure de- 
scribed in Agafonovaet al. (2005). However, the spectral shape 
can be probed by the 'trial and error' method. For example, the 
spectrum labeled as SI in Fig. 2 delivers a better fit to the ob- 
served line intensities compared to HM96. The synthetic line 
profiles for the subsystems A and B calculated with physical 
parameters obtained with spectrum SI (Col. 3 and 5 in Table 2) 
are shown by the smooth lines in Fig. 2. Part of the C IV ab- 
sorption at —30 km s~^ < u < 15 km s~^ was not included in 
the calculations because of possible input from the subsystem 
C. 

As for the highly ionized subsystem C, its hydrogen lines 
are hidden in the Hi absorption lines of subsystem B. The 
upper Umit for the Hi column density along with the col- 
umn densities for the Cm 977, Civ 1548, 1550, Oiv 832 
and Ovi 1031, 1037 lines identified in this subsystem (in- 
tensities at the expected positions of Cm and Oiv may be 
caused by unidentified blends) are given in C0I.6 of Table 2. 
The FWHM of the O vi lines is 24 km s~^ which corresponds 
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to Tkin ~ 1 .6 X 10^ K if the line broadening is entirely thermal. 
According to Sutherland & Dopita (1993), maximal output of 
the O VI ion in the case of collisional ionization is reached at 
Tkin = 3 X 10^ K and decreases quickly at lower tempera- 
tures. Thus, this subsystem may be produced either by gas that 
has been (shock)-heated and is now cooling and recombining 
(non-equilibrial ionization) or by low-density gas in photoion- 
ization equilibrium with an ionizing background hard enough 
to account for the observed coluimi density of O vi (e.g. such 
as HM96orSl). 



In general, the parameters obtained for the whole Zg,hs = 
1.8073 absorption complex - high metalUcity and low ioniza- 
tion compact clouds (L < 0.2 kpc) embedded in hot highly ion- 
ized gas - resemble gas observed in a bursted out superbubble 
by Heckman et al. (2001). The suggested vicinity to a galaxy 
makes it interesting to test an ionizing background with signifi- 
cant input from stellar radiation, i.e. a spectrum which is softer 
at > 3 Ryd compared to the QSO/AGN-dominated spectrum 
of HM96. The calculations show that any type of softer spectra 
can unambiguously be ruled out, since they poorly reproduce 
the observed line intensities and, moreover, lead to a relative 
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overabundance of carbon compared to silicon, ([C/Si]> 0). As 

far as we know, such an overabundance has never been mea- 
sured in an H II region and is not expected theoretically. 

4.2. Absorption system atz^^s = 1-7301 

The system exhibits metal lines both in low (Si III 1206, O III 
832) and high (Civ 1548, 1550; Ovi 1032, 1037) ionization 
stages (Fig. 3). The parameters obtained from the MCI calcu- 
lations with the HM96 spectrum are given in Col.2 of Table 3. 
With this ionizing background, the observed ratio of O iii 832 



and Ovi 1032, 1037 lines is reproduced only marginally: a 
good fit to the O III line is accompanied by the underestimation 
of O VI (Fig. 3, dashed lines). Unfortunately, low S/N at the 
positions of these lines does not allow us to unambiguously re- 
ject the adopted spectrum. However, there is another inconsis- 
tency - the extremely large linear size of the absorber, L ^ 200 
kpc. Linear sizes of this order of magnitude (hundreds of kpc) 
are expected for the filamentary large-scale structure, but the 
gas in the 2;abs= 1.7301 system is too dense (number density 
3 X 10"'' cm"'^ , overdensity (5h ~ 80) to be attributed to a fila- 
ment. According to both observations (Chen et al. 1998, 2001) 
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Table 3. Physical parameters of the ^abs = 1.7301 and 1.6131 metal absorbers towards HS 0747+4259 derived by the MCI procedure with the 
Haardt & Madau and modified UV background spectra (marked, respectively, by HM96 and SI in Fig. 2). Column densities are given in cm~^ 



2abs =1.7301 Zabs =1.6131 



Parameter" 


HM96 


SI 


HM96 


SI 


(1) 


(2) 


(3) 


(4) 


(5) 


TT 

Uo 


1 1 


3. /hL— z 


> 2.8E-1 


1 ITU 1 

1.3c,— 1 


AT 

/Vh 


1 .ynzu 


ft CC 1 ft 


c nu 1 ft 


O CTI? 1 ft 

Z.Dliiy 


(Jv, km s~^ 


24. y 


24.5 


26.4 


24.5 




U. /z 


U.Do 


U. 


U.04 




o.lil— / 


1 AC A 

l.Oil— 


^ 1 C A 

o.lr:,— 










<,z.uc.— o 






Z.oJi— O 


D.Oc-— O 


i.oJj/— J 




^Si 


1 ATh 7 


Z. / Cj— 1 






[7^1 


— z.ou 


-2.18 


1 on 


1 fin 








< — i.O 


< — i.i 


[7^1 


—9 9A 


— 1 Q4 




— 1 17 

— 1. -L i 


f7o.l 


9 ^7 
— Z.o / 


9 11 






/VCH ^\ 


1 .yc i u 


I i.y ziz u.4:)ciu 


7 M 


l/.Z HZ Z.UIIiiH- 


AT(f^ tttA 
iV [y. lllj 


A fid ^ 


\o.x) ± U.o jc-l J 




l^o.o ±1.0 JiiiZ 


i V \\J III ) 


1 y1 
i.Jr!,14 


yL.o ± U.4jrli4 


-5 qt::i oft 




A^CSilll) 


1.7E12 


(1.7±0.2)E12 






A^(C IV) 


3.1E13 


(3.10±0.15)E13 
2.4E14'' 


1.2E13 


(l.lib0.1)E13 


NiO IV) 


2.3E14'' 






Ar(SiiV) 


9.0E11 


(8.3±3.0)E11 






iV(Nv) 






<1.6E13 


<1.6E13 


N{0 VI) 


2.8E13 


(4.0± 1.0)E13 


1.3E14 


(1.30ib0.25)E14 


(T>,K 


4.0E4 


3.9E4 


> 4.7E4 


> 4.5E4 


ran, cm ^ 


3.2E^ 


5.6E^ 


< l.OE^ 


< 2.4E-^ 


L, kpc 


200 


56 


> 190 


^ 35 



"Zx = A^x/A^h; [.^x] = log(Arx/A^H) - log(iVx/A^H)o; Calculated using the 
corresponding velocity and density distributions estimated from the metal profiles 
marked by the horizontal bold lines in Figs. 3 and 4 



and theoretical predictions (Dave et al. 1999), strong Lya ab- 
sorbers [A'^(Hi) > 10^^ cm~^ ] arise at impact parameters to a 
galaxy of r < 50/i^^ kpc. Thus, the linear size of the Zgx^s = 
1.7301 absorber may be similar. 

The ionizing spectrum of HM96 was modified in a way 
to enhance the O Vi/O iii ratio and to increase the gas num- 
ber density [i.e. to decrease the mean ionization parameter Uq, 
see eq.(l)]. The fractional abundance of Oiii decreases with 
the enhancement of the intensity at 3 < < 4 Ryd and the 
fraction of O VI becomes larger with higher intensity at > 4 
Ryd. One of the acceptable spectral shapes is shown in Fig. 2 
by the short-dashed line labeled as SI. Results obtained with 
SI are given in Col.3 of Table 3, and the synthetic profiles are 
plotted in Fig. 3 by the smooth lines. Note that with the same 
normalizing intensity J912 the spectrum SI gives nearly four 
times smaller linear size than the spectrum of HM96. 

It should be stressed that the shape SI is only one of the 
possible solutions for the z^'os = 1.7301 system since the low 
S/N of most of the spectral data allows us to vary the intensity 
at £ > 3 Ryd over a broad range. Thus, the exact spectral 
shape cannot be recovered. Nevertheless, we can conclude that 
the absorption fines observed in this system favor spectra with 
increased intensity at £ > 3 Ryd as compared to the HM96 
model. 



4.3. Absorption system atz^hs = 1-6131 

This highly ionized system exhibits several hydrogen lines 
along with a weak C IV doublet and strong lines of O VI 103 1, 
1037 (Fig. 4). At the expected position of the Siiv 1393 fine 
a clear continuum is present. The observed intensity at the po- 
sition of the C III 977 fine can be used only as an upper fimit 
on the C III absorption since this line lies in a noisy part of the 
spectrum and is partly contaminated by LyS from the Zahs = 
1.7219 system. This upper limit, along witfi a safe upper fimit 
on the abundance ratio [0/CJ < 0.5 known from measurements 
in Galactic and extragalactic H II regions and in metal poor halo 
stars (Henry et al. 2000; Akerman et al. 2004), are utilized to 
estimate the physical parameters in the absorber. The parame- 
ters obtained from the MCI calculations with the HM96 spec- 
trum are listed in Col.4 of Table 3, whereas those obtained with 
the spectrum SI (Fig. 2) are given in Col.5. The synthetic line 
profiles are identical in both cases and are shown by the smooth 
lines in Fig. 4. 

The spectrum of HM96 leads to a linear size of hundreds of 
kpc and, hence, places the absorber in some filamentary struc- 
ture. The spectrum SI is harder at > 3 Ryd and provides 
a fine-of-sight size several times smaller, which allows us to 
attribute the absorber to a galactic halo. 
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Fig. 5. Same as Fig. 1 but for the «abs = 1.5955 system. Smooth curves are the synthetic spectra calculated with the ionizing background S2 
(Fig. 2). The zero radial velocity is fixed at z = 1.5955. The physical parameters corresponding to the background S2 are listed in Table 4, 
Col.3, whereas those obtained with HM96 are in Table 4, Col 2. 



Both results are obtained assuming photoionization equi- 
librium for all observed ions. However, collisional ionization 
of O VI should be tested as well. From the apparent FWHMs 
of the available lines: H i (60 km s~^ ), C iv (25 km s~^ ), and 
O VI (35 km s~^ ), and assuming pure thermal line broaden- 
ing, we obtain kinetic temperatures of Thi = 7.6 x 10^ K, 
Tciv = 1.4 X 10^ K, and Tqvi = 3.8 x 10^ K. Such diverse 
temperatures carmot occur in the same gas and thus turbulent 
broadening must be significant. The possibility of coUisionally 
ionized O VI (with corresponding broad and shallow hydrogen 
lines hidden in the observed profile of H i) surrounds cooler gas 
producing H I and C iv absorption seems to be quite unproba- 



ble since the lines of C iv and O vi are well aligned in velocity 
space and have similar profiles. The case of non-equilibrium 
ionization can be ruled out as well since the recombination time 
of O VI is much shorter than that of hydrogen: for example, at 
Tkin ~ 3 X 10^ K it is about 60 times faster (Osterbrock 1974). 
Then the bulk of O vi should come to photoionization equilib- 
rium prior to H I and we would observe a weak and shallow H I 
Lya along with strong O VI fines. This contradicts the observed 
profiles in the 2;abs =1.6131 system. Therefore photoionization 
equilibrium is the most probable assumption. 

It should be noted that O VI lines trace more rarefied gas not 
seen in C IV and, hence, can be broader. On the other hand, the 
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Fig. 6. Same as Fig. 1 but for the Zabs = 1.5401 system. The zero radial velocity is fixed at 2: = 1.5401. The smooth curves are the synthetic 
spectra calculated with the ionizing background S2 (Fig. 2). The physical parameters corresponding to the background S2 are listed in Table 4, 
C0I.5, whereas those obtained with HM96 are in Table 4, Col.4. The central positions of blends are indicated by tick marks. 



H I profile can be wider than O VI since some low density re- 
gions produce feeble O Vl absorption not distinguishable from 
the noise. Thus, we can estimate the lower limit on the tur- 
bulent velocity dispersion from the comparison of the H i and 
O VI lines: aturb > 19 km s~^ . 

Taking into account such high turbulent velocities, the 
placement of the Zabs = 1-6101 absorber in a galactic halo, i.e. 
within boundaries where it remains gravitationally bound to a 
galaxy, seems to be more reasonable than its origin in the fila- 
mentary structure with no significant sources of the line broad- 
ening apart from the Hubble expansion. Thus, ionizing spectra 
that are harder atE > 4 Ryd as compared to that of HM96 are 
preferable for this absorber. 

A few words should be said about the high upper limit 
on the nitrogen abundance since measurements usually show 
[N/C]< -0.5 for metal poor gas (Cenfiirion et al. 2003). 
Unfortunately, the N v 1238 line in the absorber under study 
lies in a very noisy part of the Keck spectrum with many broad 
continuum undulations of unidentified nature. It is possible that 
this low contrast feature is not completely due to N V 1238 ab- 
sorption. 



4.4. Absorption system Zabs = 1-5955 

Absorption lines identified in this system are shown in Fig. 5. 
Note the very strong lines of O Vl 1032, 1037 observed together 
with lines of Si ill 1206 and Si iv 1393 (Si IV 1402 is blended). 
The physical parameters derived by the MCI assuming the 
HM96 ionizing background are listed in Col. 2 of Table 4, and 
the corresponding synthetic profiles are shown by dashed lines 
in Fig. 5. This ionizing spectrum is not optimal: for the distri- 
bution of ionization parameters determined by the C III, C IV 
and O VI lines it significantly underproduces intensities of siU- 
con lines (or, equivalently, delivers the ratio Si/O significantly 
higher than solar, which is inconsistent with measurements in 
Hii regions) and it gives too large of a linear size. For this 
system non-equilibrium ionization can also be ruled out: the 
broadening of the Cm 977, Civ 1548, 1550, Ovi 1031, 1037 
and hydrogen lines is similar. This points to significant contri- 
butions from turbulent velocities and, hence, to temperatures 
of about 10^ K (see previous subsection). Moreover, there is a 
safe upper limit on the nitrogen abundance [N/C]< —0.5 ob- 
tained from the Nv 1238 line: since Nv ttaces the same gas 
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as O VI, any 'overionization' of this gas would result in an arti- 
ficial enhancement of the nitrogen abundance calculated under 
the assumption of equifibrium ionization. Thus, we look for a 
more appropriate ionizing background. 

Taking into account information obtained with the HM96 
spectrum and using the technique described in Reimers et al. 
(2005) and Agafonova et al. (2005), we searched for a back- 
ground that would ensure a higher ratio of fractional ioniza- 
tions Tsiiv/Tciv at the ionization parameter U determined 
from the C III, C IV and O VI absorption. Additionally, this U 
should be low enough to produce a physically reasonable lin- 
ear size for the absorber under study. Other constraints to the 
solution were the inequalities [Si,0/C]< 0.5. 

The resulting spectrum is shown in Fig. 2 by the dashed line 
and marked as S2. It is significantly harder in the whole range 
E > 1 Ryd compared to both the initial spectrum of HM96 
and the spectrum SI used for the absorbers described in the 
preceding subsections. The physical parameters obtained from 
the MCI calculations with the S2 spectrum are listed in Col. 3 of 
Table 4, and the corresponding synthetic profiles are plotted by 
the smooth fines in Fig. 5. Note that the derived element abun- 
dance ratios [Si,0,N/C] are consistent with measurements in 
Galactic and extragalactic H ii regions (e.g., Matteucci 2003). 

The assortment and column densities of lines observed in 
the 2abs= 1.5955 system resemble those identified in some 
highly ionized high- velocity clouds (HVCs) at z = (CoUins 



et al. 2004). Perhaps we are observing in the Zabs = 1.5955 ab- 
sorber a high-redshift analog of the local HVCs. 

One important difference between the spectra S2 and SI 
should be emphasized. The shape S 1 represents only one pos- 
sible solution from a wide range of acceptable shapes since the 
low S/N for lines observed in the .2abs= 1.8073, 1.7301 and 
1.6131 systems does not permit us to recover the SED with 
reasonable accuracy. Conversely, lines observed in the z^hs = 
1.5955 system can be self-consistently described only with a 
quite narrow range of possible ionizing spectra and, hence, the 
spectrum S2 is determined much more accurately. In this con- 
text it is interesting to test the shape S2 on the absorption sys- 
tems at Zabs = 1.8073 and 1.7301 exhibiting many metal lines. 
The result is negative - this spectral shape is inconsistent with 
line strengths observed in these systems. Thus, the ionizing ra- 
diation at 2 < 2 seems to undergo strong spectral variations. 

4.5. Absorption system a? Zabs = 1.5401 

This absorption system is seen in H I Ly/3 and Ly7 and in metal 
lines shown in Fig. 6 (Lya and Si iii 1206 fall in the wavelength 
gap between the HST and Keck spectra). Note the unusually 
narrow C iv lines with FWHM = 12.5 km s~^ . 

The parameters derived from the MCI calculations with the 
HM96 ionizing spectrum are listed in Col.4 of Table 4. The 
synthetic profiles are shown in Fig. 6 by dashed lines. This 
spectrum is not completely consistent with the observed line in- 
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Table 4. Physical parameters of the Zabs = 1.5955, 1.5401 and 1.4649 metal absorbers towards HS 0747+4259 derived by the MCI procedure 
with the Haardt & Madau and modified UV background spectra (marked, respectively, by HM96, SI, and S2 in Fig. 2). Column densities are 
given in cm~^ 
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"Zx = Ny^/Nu; [Zy.] = log{N:si/NH) - \og{Nx/NH)Q; ''calculated using the velocity and density distributions 
estimated from the metal profiles marked by the horizontal bold lines in Fig. 5 (Zabs =1.59) and in Fig. 6 (2abs =1.54) 



tensities: it overestimates C III (the intensity of C III should be 
lower than the apparent intensity due to blending with Lyd from 
the Zahs = 1.6131 system) and significantly underestimates the 
profiles of C IV 1548, 1550 in the center. 

The spectra S 1 and S2 were tried as well. S 1 was rejected 
for the same reason as the spectrum of HM96 - poor reproduc- 
tion of carbon lines. The spectrum S2 turned out to be consis- 
tent with the data. The physical parameters obtained with this 
background are listed in Col.5 of Table 4, and the correspond- 
ing synthetic line profiles are plotted in Fig. 6 (smooth lines). 

This absorption system has characteristics which are quite 
unusual for intergalactic absorbers: metallicity of 0.6 solar and 
a sub-kiloparsec linear size. Probably it was ejected from a 
nearby galaxy located transverse to the line of sight. Note that 
the hard ionizing continuum derived for both z^hs = 1 .5955 and 
-^abs = 1.5401 absorbers (separated by 6500 km s~^ ) may be 
considered as a hint to the presence of an AGN in the vicinity of 
these systems (cf. S2 and the AGN-type spectrum of Mathews 
& Ferland in Fig. 2). 

4.6. Absorption system atzgi^s = 1-4649 

This system is almost identical to that at 2:abs= 1.6131 (the 
lines Nv 1238, 1242 fall in the gap between Keck and STIS 
spectra). The physical parameters obtained from calculations 



with the HM96, SI and S2 spectra are listed in Col.6, 7 and 8 
of Table 4, respectively. The synthetic profiles coincide for all 
three cases and are shown by the smooth lines in Fig. 7. The 
system exhibits a similar high turbulent velocity as found at 
-^abs = 1.6131. Therefore its origin in a filamentary structure is 
less favorable than an origin in a galactic halo and, hence, the 
harder spectra SI and S2 are preferable to that of HM96. 

4.7. Frequency distribution of Oy\ absorbers 

The detection limit of O VI A 1302 lines in the spectrum of 
HS 0747+4259 corresponds to iV(Ovi) > 2 x lO^^ cm^^ 
(W^rest ^ 30 mA). In addition to the 6 absorbers described 
above, a further 10 metal systems in the redshift range 1.07 - 
1.87 exhibit O vi A1031 absorption with the second line A1037 
blended with Lya forest absorption. However, the identifica- 
tion of O VI in these systems is certain due to the simultaneous 
detection of the C iv and O iv A787 lines (in absorbers with 
-^abs > 1.7). Table 5 lists the redshifts and column densities for 
the detected O IV lines. 

With the absorption distance of AX = 1.27 (go = l/2)\ 
the number of encountered O vi absorbers translates into an 
absorber frequency of AAovi/AX = 13. This is almost two 
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Table 5. O vi column densities 





^ahs 


iV(0 VI), cm-2 


Comments 


1 


1.8617 


(3.2 ± 0.8)E13 


A1037 blended 


2 


1.8533 


(3.0±0.8)E13 


A1037 blended 


3 


1.8073 


(1.1 ±0.3)E14 




4 


1.7790 


(2.5±0.8)E13 


A1037 below detection limit 


5 


1.7744 


(4.0 ± 1.0)E13 


A1037 blended 


6 


1.7302 


(4.0 ± 1.0)E13 




7 


1.7154 


(2.5±0.8)E13 


A1037 blended 


8 


1.6331 


(7.0 ± 1.5)E13 


A1037 blended 


9 


1.6134 


(1.3 ± 0.3)E14 




10 


1.5955 


(4.2 ± 0.8)E14 




11 


1.5401 


(1.9±0.4)E13 




12 


1.4856 


(1.1±0.3)E14 


A1037 blended 


13 


1.4648 


(1.7±0.4)E14 




14 


1.2910 


(5.0±3.0)E13 


noisy spectrum, S/N ^ 3 


15 


1.1896 


(7.0±3.0)E13 


A1037 blended 


16 


1.0778 


(5.0±3.0)E13 


noisy spectrum, S/N ~ 3 



times higher than AJ\fovi/AX = 7 measured for the range 
1.21 < 0abs< 1.67 (AX = 0.72) in the spectrum of HE 0515- 
4414 (Reimers et al. 2001). Thus, the O VI absorber statistics at 
^abs ~ 1 -5 is subject to strong variations from sightline to sight- 
line. The mass density of the O vi absorbers, J2 -^(O 'Vi)/AX, 
shows even stronger variations: (1.2 ±0.3) x 10^^ cm~^ along 
the present sightline compared to 3.0 x 10^^ cm^^ towards HE 
0515-4414. Note, however, that the bulk of the integrated col- 
umn densities in HS 0747+4259 comes from the strong O VI 
absorbers clustered at 1.46 < Zabs < 1-63 where hard ionizing 
background was recovered. For comparison, at low redshift, 
z < 0.15, the O VI absorber frequency (VFics > 30 mA) av- 
eraged over 31 FUSE sightlines is AA/bvi/AX = 17 ± 3 
(Danforth & Shull 2005). We recognize, however, such a com- 
parison neglects the fact that the Ovi systems at z ^ 1.5 are 
mainly photoionized, while those in the local universe may 
originate in the colUsionaUy ionized intergalactic medium. 

5. Conclusions 

In this work, we investigate the O Vi absorption Une systems 
detected in the HST/STIS and Keck/HIRES spectra of the QSO 
HS 0747+4259. Six systems O VI with Zabs = 1 -46- 1 .8 1 reveal a 
sufficient number of lines to perform a quantitative analysis of 
the shape of the UV background ionizing radiation. The results 
obtained lead to the following conclusions. 

1. Spectral energy distributions with a significantly (> 2.5 
times) higher intensity at E > 4 Ryd as compared to the 
model spectrum by HM96 are preferable for the absorption 
systems studied. Similar conclusion have been obtained in 
Agafonova et al. (2005) for the absorber at z^hs = 1-9426 
towards J 2233-606. Such SEDs are suggested by the de- 
ficiency of strong Lya absorbers [AT (Hi) > 10^^ cm^^ ] 
at ^; < 2. On the other hand, the metagalactic SEDs ob- 
tained at z ~ 3 (Agafonova et al. 2005) show an intensity at 
E > 4 Ryd consistent with that of HM96. This confirms the 
conclusion of Kim et al. (2001) that absorbers with strong 
Lya lines evolve faster than weak ones. 



2. The spectral shapes recovered for the absorbers at Zabs = 
1.8073, 1.7301 and at Zabs= 1.5955, 1.5401 are different 
with the latter being significantly harder in the whole range 
E > 1 Ryd. There are theoretical predictions (Dave et al. 
1999) that in general the evolution of the Lya absorbers 
flattens at z < 1.7. However, this flattening would not lead 
to the increased ionizing flux. Thus, the observed variation 
of the spectral shape is probably of local nature and may 
be caused by the presence of a QSO/AGN close the Une of 
sight (cf. Jakobsen et al. 2003). 

3. The spectral shapes at z < 2 do not show signs of a 
significant contribution from galactic (stellar) radiation. 
According to Bianchi et al. (2001), an escape fraction 
/esc = 0.05 results in a galactic contribution to the UV 
background at 912 A equal to that of QSOs. Such galac- 
tic input would noticeably soften the ionizing spectrum in 
the range E > 3 Ryd (see, e.g., Giroux & Shull 1997). 
Taking into account the hardness of the recovered spectra 
at ii^ > 3 Ryd , the escape fraction of galactic photons 
at 2; < 2 should be well below 0.05. This is in line with 
/esc < 0.04 estimated by Fernandez-Sotoetal. (2003) from 
measurements of 27 galaxies with 1.9 < z < 3.5. 

4. In the systems considered, the main ionization mechanism 
for all observed ions including O VI is photoionization by 
the background UV radiation. The contribution of turbu- 
lence to the line broadening is significant and the tempera- 
ttires estimated from the widths (FWHMs) of the O VI lines 
are well below the range where collisional ionization is 
dominant. High O VI column densities (N > lO^"' cm^^ ) 
can originate in non-equilibrium gas which has been shock- 
heated and is now cooling and recombining (Heckman et 
al. 2002). Such conditions can be valid for very particular 
absorption system at 2;abs= 1.8073 (subsystem C), but no 
other systems show signs of being non-equilibrial. Similar 
conclusions about the ionization of high redshift O VI are 
reported by Carswell et al. (2002), Bergeron et al. (2002), 
and Levshakov et al. (2003). 
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